First-principles calculations of strontium on Si(OOl) 
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This paper reports state-of-the-art electronic structure calculations on the deposition of strontium 
on the technologically relevant, (001) orientated silicon surface. We identified the surface reconstruc- 
tions from zero to four thirds monolayers and relate them to experimentally reported data. A phase 
diagram is proposed. We predict phases at 1/6, 1/4, 1/2, 2/3 and 1 monolayers. Our results are 
expected to provide valuable information in order to understand heteroepitaxial growth of a promi- 
nent class of high-K oxides around SrTiOs. The insight obtained for strontium is expected to be 
transferable to other alkaline earth metals. 

PACS numbers: 68.43.Fg, 68.47.Fg, 71.15.Mb, 73.20.-r 



I. INTRODUCTION 

Device scaling has been the engine driving the micro- 
electronics revolution as predicted by Moore's lawi By 
reducing the size of transistors, processors become faster 
and more power efficient at an exponential rate. Cur- 
rently the main challenge in device scaling is the integra- 
tion of high-K oxides as gate oxides into silicon technol- 
ogy 

The gate oxide is the dielectric of a capacitor, which 
is used to attract charge carriers into the channel region. 
Thus a current can flow from source to drain, provided 
a voltage is applied to the gate electrode. With a thick- 
ness of only 1-2 nm^ the gate dielectric is the smallest 
structure of a transistor. As the thickness of the gate 
oxide is further reduced, its insulating property is lost 
due to direct tunneling through the ultrathin oxide. The 
results are intolerable leakage currents and a large power 
consumption. 

A remedy to this problem is the replacement of the 
current Si02 based gate oxides with an insulator having 
a larger dielectric constant, a so-called high-K oxide. A 
high-K oxide gate with the same capacity as an ultrathin 
Si02 based one will be thicker and should therefore ex- 
hibit smaller leakage currents due to direct tunneling. 
The integration of new oxides into the semiconductor 
technology has, however, proven to be a major problem. 
Hence an enormous research effort is underway to under- 
stand growth of high-K oxides onto silicon. 

Currently, Hf02 and Zr02 are the main contenders for 
the first generation of high-K oxides to be introduced in 
fabrication. These oxides still exhibit an interfacial Si02 
layer and therefore do not form a direct interface with sil- 
icon. As sealing proceeds, an interfacial Si02 layer can- 
not be tolerated anymore. The existence of an atomically 
abrupt interface between silicon and a high-K oxide has 
been demonstrated by McKee et alii^ for Baa;Sri_a;Ti03 
on Si(OOl), after an epitaxial relationship has been re- 
ported in the late 80's.^ 

A detailed understanding of metal adsorption is crucial 



to control oxide growth on Si. The growth process is 
guided by the sequence of structures that develop as the 
metal is deposited on the silicon surfacei^ The nature of 
these structures as well as the interface between Si and 
a high-K oxide is, however, still under debate. 

The adsorption of the alkaline earth metals Sr and Ba 
on Si(OOl) has been extensively studied. Most of the 
studies of Sr on Si(OOl) are diffraction studies such as low- 
energy electron diffraction (LEED)i'2*2iifliiiiiS and reflec- 
tion high-energy electron diffraction (RHEED)iiili^ 
or scanning tunneling microscopy (STM)^*^*!^ experi- 
ments. The STM studies have been most valuable be- 
cause they contribute detailed real-space information on 
the atomic scale. Simflar LEED,iLi2iia RHEEDSa and 
grpjy^2ij22j23j24j25 ^^^1 as X-ray photoemission stud- 
ie^iSi^ have been performed for Ba. X-ray standing wave 
experiments provide valuable restrictions on the struc- 
tures with coverages of 1/2 and 1/3 monolayer (ML)i21 
The photoemmision (XPS)iiiS£ studies provide insight 
into the ionization state of Sr, and show a qualita- 
tive change of the Fermi-level pinning as a function of 
coverage 
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Diffraction studies suffer from the fact that they av- 
erage over several structures and terraces. Here STM 
experiments provide valuable clues. One of the major 
experimental difficulties is the determination of the cov- 
erage at which the data are coUcctediii 

Theoretical investigations of isolated Ba atoms ad- 
sorbed on Si(OOl) have been performed by Wang et ali^S 

In this work we address the deposition of Sr on Si(OOl) 
using state-of-the-art electronic structure calculations. 
We attempt to provide a complete set of adsorption 
structures, their energetics, chemical binding and elec- 
tronic structure. We will categorize the reconstructions 
by pointing out the driving forces that lead to the various 
ordered structures. This provides a unified picture of Sr 
adsorption from low coverage up to 4/3 monolayers. 
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II. COMPUTATIONAL DETAILS 

The calculations are based on density functional 
theorji^Si^ using a gradient corrected functional The 
electronic structure problem was solved with the projec- 
tor augmented wave (PAW) method;^ an all-electron 
electronic structure method using a basis set of plane 
waves augmented with partial waves that incorporate the 
correct nodal structure. The frozen core states were im- 
ported from the isolated atom. For the silicon atoms we 
used a set with two projector functions per angular mo- 
mentum for s and p-character and one projector per an- 
gular momentum with o?-character. The hydrogen atoms 
of the back surface had only one s-type projector func- 
tion. For strontium we treated the 4s and 4p core shells 
as valence electrons. Per angular momentum we used 
three s-type and two p- and d-type projector functions. 
The augmentation charge density has been expanded in 
spherical harmonics up to = 2. The kinetic energy cut- 
off for the plane wave part of the wave functions was set 
to 30 Ry and that for the electron density to 60 Ry. 

A slab of five silicon layers was used as silicon sub- 
strate. Wang et aliS report that the adsorption energy 
of a Sr atom on the surface changes by 0.05 eV when be- 
tween a 4-layer slab and a 6-layer slab of silicon. In our 
calculations the energy per additional silicon atom agrees 
to within 0.06 eV with that of bulk silicon between a 4- 
layer and a 5-laycr slab. The dangling bonds of the unre- 
constructed back surface of the slab have been saturated 
by hydrogen atoms. The lateral lattice constant was cho- 
sen as the experimental lattice constant a = 5.4307 A of 
siliconi^ The slabs repeat every 16 A perpendicular to 
the surface, which results in a vacuum region of 9.5 A for 
the clean silicon surface. 

The Car-ParrincUo ab-initio molecular dynamics^ 
scheme with damped motion was used to optimize the 
electronic and atomic structures. All structures were 
fully relaxed without symmetry constraints. The atomic 
positions of the backplane of the slab and the terminating 
hydrogen atoms were frozen. 

Many of the Sr adsorption structures are metallic 
which requires a sufficiently fine grid in k-space. We used 
an equivalent to eight by eight points per (1x1) surface 
unit cell. This value has been chosen after careful conver- 
gence tests for surface structures, bulk silicon and bulk 
Sr silicides (Fig.^). In cases where this k-mesh is incom- 
mensurate with the size of the unit cell we selected the 
closest, finer commensurate k-mesh. 

For metallic systems, the orbital occupations were de- 
termined using the Mermin functionalSS, which produces 
a Fermi-distribution for the electrons in its ground state. 
The electron temperature was set to 1000 K. In our case 
this temperature should not be considered as a physi- 
cal temperature but rather as a broadening scheme for 
the states obtained with a discrete set of k-points. The 
Mermin functional adds an entropic term to the total 
energy, which is approximately canceled by taking the 
mean of the total energy U{T) and the Mermin- free cn- 
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FIG. 1: K-point convergence: Energy in eV relative to the 
converged value versus characteristic k-point spacing!— From 
top to bottom: bulk silicon per atom; adsorption energy per 
Sr-atom for a coverage of 1/2 ML; bulk silicides SrSi Cmcm, 
SrSi2 lAi/amd and SrSi2 P4332 per atom. Our surface calcu- 
lations used a k-point spacing of 0.2 A^^ or the closest com- 
mensurate mesh. 



Energy [H] 


So [Si] 




-4.0036 


So[Sr] 




-31.1441 


i5o[5 layer- 


Si-slab] 


-21.1140 


£"0(4 layer- 


Si slab] 


-17.1083 



TABLE I: Reference energies used in this paper without 
frozen core energy. See text for details of the calculation. 

crgy F{T) = U{T) - TS{T) as proposed by Gillan»2a 

C/(T = 0)«i(F(T) + [/(T)). (1) 

In order to express our energies in a comprehensible 
manner, we report all energies relative to a set of refer- 
ence energies. This set is defined by bulk silicon and the 
lowest energy polymorph of SrSi2 (P4332). The refer- 
ence energies are listed in Tab. The reference energy 
Ef) [Sr] for a Sr atom, corresponding to the coexistance of 
bulk silicon and bulk SrSi2, is extracted from the energy 
E[SrSi2] of the disilicide calculated with a (8 x 8 x 8) 
k-mesh and the reference energy of bulk silicon i?o [Si] as 

£;o[Sr] =£;[SrSi2] -2i;o[Si]. (2) 

The bulk calculation for silicon was performed in the 
two atom unit cell with a (10 x 10 x 10) k-mesh and at 
the experimental lattice constant of 5.4307 Ai^ 

For the surface calculation, we always subtracted the 
energy of a clean (4 x 2) silicon surface of the same slab 
thickness, to account for slab including hydrogen termi- 
nation. For the surfaces with half a monolayer coverage 



3 



of silicon we assumed that the corresponding reservoir for 
the sihcon atoms is a sihcon terrace. Hence, the reference 
energy for the 4.5 layer silicon slab is the average energy 
of a 4-layer slab and a 5-layer slab. The terrace energy 
itself does not enter the limit of an infinitely dilute step 
density. 

In some of our structures the choice of unit cell has an 
impact on the dimcr buckling. Wc estimated the energy 
of a buckling reversal from the energy difference of a (2 x 
2) and a (5 x 2) supercell. The cell with an odd number 
of dimers contains one buckling reversal. The calculated 
energy for such a buckling reversal is 0.06 eV. 

III. BULK SILICIDES 

Before studying the adsorption of Sr on silicon wc in- 
vestigated the bulk silicidcs of Sr. The energetics of bulk 
silicidcs provide us with the driving force to go from or- 
dered surfaces structures to silicidc grains on the surface. 
Our calculations on early transition metals on silicon in- 
dicate that silicide formation is a major problem for layer- 
by-layer growth of an oxide^ 

The binding characteristics of the bulk silicidcs pro- 
vide us with insight into the favored structural templates 
which might be anticipated for the Sr-covcrcd silicon sur- 
face. 

Sr-silicides are typical Zintl compounds. According 
to the Zintl-Klemm conccptjiS, atoms with an increased 
number of electrons form similar structures as atoms with 
the correspondingly increased atomic number. Conse- 
quently, a charge transfer of one electron to silicon will re- 
sult in a preferred bonding environment similar to phos- 
phorous with three or five covalent bonds. Addition of 
two electrons will result in chain-like structures like sul- 
fur. After addition of three electrons one anticipates for- 
mation of dimers and once four electrons are transfered, 
isolated ions are expected. In other words: for every 
added electron one covalent bond will be missing. 

Due to the large difference in electronegativity. Sr for- 
mally donates its two valence electrons to the silicon sub- 
strate. 

• In SrSi2 one electron is transfered per silicon atom. 
Hence three-fold coordinated silicon networks are 
formed as shown in the top two structures of Fig. |21 

• Two electrons are transfered in SrSi so that the 
silicon network is similar to that of elemental sul- 
fur with two-coordinated silicon atoms forming Si 
chains as seen in the middle left panel of Fig. |21 
The middle right panel shows another modification 
of SrSi. The average number of covalent bonds per 
silicon atom is, however, still two. 

• In SrsSia there is a charge transfer of 10 electrons to 
three silicon atoms, which can be used to form two 
Si^~ ions and one Si^~ ion. The two Si'^" combine 
to form dimers and the Si^~ is no more able to form 




FIG. 2: Bulk silicide structures. Top left: SrSi2 (/4i/amd)ii; 
Top right: SrSi2 (P4332)^; Middle left: SrSi {Cmcm)^- Mid- 
dle right: SrSi (Immm)^; Bottom left: SrsSia (74/mcm)^; 
Bottom right: Sr2Si {Pnma)~. The large, dark spheres rep- 
resent Sr atoms, the smaller, light spheres Si atoms. Energies 
are listed in Tab.|lTl 

covalent bonds. Hence we observe an equal number 
of Si dimers and single Si ions in the structure of 
SrsSis as shown in the lower right panel of Fig.|21 

• In Sr2Si four electrons are transfered to each sili- 
con atom. As a consequence, the Si atoms in the 
stucture on the bottom right of Fig. |21do not form 
covalent bonds. 

We find SrSi2 (P4332) to be the most stable phase 
of silicidcs per Sr atom (Tab. EJ. Therefore we have 
chosen this material to define, together with bulk Si, the 
reference energy for Sr. 

As a side remark, we note that a lower energy of a Sr 
atom in a bulk silicide compared to the adsorbed Sr on 
the surface does not automatically indicate the forma- 
tion of silicide grains during growth: The silicide forma- 
tion may be suppressed by the strain due to an epitax- 
ial constraint by the silicon lattice constant. Thus the 
formation of silicidcs is expected to be delayed for thin 
films, because the bulk silicidcs have a large mismatch 
with the silicon substrate. This argument does not refer 
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E[Sr] [cV] 



SrSi2 (P4332) 0.00 
SrSi2 {I4i/amd) 0.01 
SrSi [Cmcm) 0.09 
SrSi (Immm) 0.20 
SrsSis (/4/mcm) 0.39 
SraSi {Pnma) 0.45 



TABLE II: Energies per Sr atom of bulk silicides relative to 
our reference energies. 

to the thermodynamic equilibrium of large samples, but 
it indicates that nuclcation of silicidc grains will have to 
overcome a large barrier. 



IV. THE CLEAN SILICON SURFACE 

The clean (001) silicon surface has a c(4 x 2) dimer-row 
reconstruction. We briefly summarize the driving forces 
towards this reconstruction in order to understand the 
adsorption structures of Sr on silicon. 

An unreconstructed surface of silicon (001) is termi- 
nated by a square, (1x1) array of atoms. Each silicon 
atom on the surface is connected by two bonds to the 
subsurface. Consequently, there are two half-occupied 
dangling bonds on each silicon sticking out of the sur- 
face. 

Each pair of surface silicon atoms forms a dimcr bond 
which saturates one of the two dangling bonds on each 
atom. The dimers arrange in rows. This is the so-called 
(2x1) dimer row reconstruction which results in an en- 
ergy gain of 0.65 eV per (1 x 1) unit cell according to our 
calculations. 

In a second reconstruction, both electrons in the dan- 
gling bonds localize on one atom of each dimcr, result- 
ing in a dimer buckling. The buckling is driven by 
the fact that a sp'^ hybridization is favored for a five- 
electron species such as the negative Si atom, while a 
sp^ hybridization is favored for a 3-electron species such 
as the positive silicon atom. The sp^ hybridization in 
turn favors a planar bonding environment, whereas three- 
coordinated sp'^ bonded atoms form an umbrella like en- 
vironment. In the buckled-dimcr row reconstruction, the 
electrons are localized on the atoms sticking out farthest 
from the surface. The dimer buckling can be considered 
as a Peierl's distortion which splits the half-filled energy 
bands resulting from the dimcr bonds into a filled and an 
empty band, with a band gap in between. The energy 
gain due to this distortion is 0.12 eV per (1 x 1) cell. 

The energies quoted here are in reasonable agree- 
ment with previously published LDA pseudo-potential 
calculations 4i The dimer reconstruction can be consid- 
ered to be fairly stable. Even at 1500 K only 3% of the 
dimer bonds are broken as estimated from the Boltzmann 
factor with AE = 0.65 + 0.12 cV. 



Two neighboring buckled dimer rows interact only 
weakly. We obtain an energy difference of 1.2 meV per 
dimer between the c(4 x 2) reconstruction with anti- 
parallel buckling and the p{2 x 2) reconstruction with 
parallel buckling on neighboring dimer rows. This in- 
dicates that the buckling patterns of different rows are 
fairly independent of each other. 

Within a row, however, the buckling of the dimers is 
coupled in an anti-correlated manner. This can, at least 
partly, be explained by the fact that the lower silicon 
atom of a dimer pushes the two adjacent subsurface sil- 
icon atoms apart. For the equivalent silicon atom of the 
next dimer, it is therefore favorable to be in the higher, 
sp^-like configuration!^ 

There has been an intense debate as to whether there 
is dimer buckling or not. STM images reveal a 2 x 1 struc- 
ture. They exhibit the buckling only at rather low tem- 
peratures and near defects. The theoretical predictions 
depend strongly on the approach chosen (cluster calcula- 
tions with configuration interaction or density functional 
calculations with periodic boundary conditions). The 
most conclusive results have been produced by Quan- 
tum Monte Carlo (QMC) simulations^ indicating that 
the buckling is present and density functional calcula- 
tions just overestimate the energy difference. 

The fact that STM experiments cannot resolve the 
dimer buckling may be due to thermal averaging of the 
two buckled configurations. We believe that the mecha- 
nism is due to the migration of a soliton-Iike defect in the 
anti-correlated buckling pattern of a dimcr row. The cal- 
culated energy for this defect is 0.06 eV (see section 111)1 . 
Thus we predict a concentration of one such defect per 
11 dimers at room temperature. Typical tunnel currents 
in STM experiments are around 1 nA, which corresponds 
to six electrons per nanosecond. 

Thus even a soliton migration barrier as large as 
0.05 cV would imply that the buckling changes once dur- 
ing the transfer of a single electron. These estimates 
should be taken with caution, since the small energy dif- 
ference of 0.06 eV carries a large relative error bar. 

From the comparison between DFT and Quantum 
Monte Carlo calculations one can deduce an error bar of 
0.05 eV per (1 x 1) unit cell due to electron corrclationsi^ 
We verified that this correction does not qualitatively af- 
fect the findings reported in this paper. 



V. ISOLATED SR ON SILICON 

In order to determine the low-coverage limit of Sr- 
adsorption wc investigated the energy as a function of 
the lateral position of a Sr ad-atom on the surface. The 
total energy as a function of the lateral position of Sr 
is obtained by constraining the lateral movement of the 
Sr-atom relative to the rigid back plane of the slab. The 
calculations were performed in a (4 x 4) surface super cell. 
We considered the high symmetry points of the (2x1) 
surface shown in Fig. |2| and the mid-point between the 
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FIG. 3: Top view of the Si(OOl) surface and the four high 
symmetry positions spanning the surface irreducible (2x1) 
unit cell. The dimer buckling is not shown. The energies are 
hsted in Tab. ITTTl 



Position 



Ag [eV] 



A 
B 
C 
D 

A-D midpoint 



0.00 
0.55 
0.75 
0.29 
0.61 



TABLE III: Relative energies of isolated Sr on the Si(OOl) 
surface at the high-symmetry points of the surface irreducible 
unit cell as well as the A-D midpoint. The labels refer to 
FigH 



local minima. The energies of the high symmetry posi- 
tions are given in Table UTTl 

Sr has the global minimum at position A as defined 
in Fig. 13 Sr is located in the trench between the dimer 
rows and in the center of 4 surrounding dimers. The Sr 
atoms are slightly elevated above the plane of the surface 
dimers. 

A metastable position, D of Fig. O is located in be- 
tween two dimers on top of a dimer row. It is 0.29 eV 
higher in energy than the global minimum. We will see 
that the structures A and D are repeating motifs in a 
range of different adsorption structures. 

The diffusion of Sr on the silicon surface proceeds 
about equally fast parallel and perpendicular to the 
dimer rows, with a slight preference for the parallel direc- 
tion. The diffusion barrier along the valley is equal to the 
energy difference between sites A and B, namely 0.55 eV, 
the one across the row is 0.61 eV and is estimated by the 
midpoint between the sites A and D 

It should be noted that in our analysis we ignored the 
reduced symmetry due to dimer buckling. As a result, 
different versions of the high symmetry points quoted 
here exist with slightly different energies. For the struc- 
ture A we found two versions which differ in energy by 
0.15 eV. In these cases the lowest energy structure has 
been chosen. 

In contrast to Sr, there has been a lot of work related to 
isolated Ba atoms adsorbed on Si(001)i2ii2S Our results 
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FIG. 4: Schematic representation of the isolated Sr ad-atom 
at position A. The filled circle represents the Sr ad-atom, 
the rectangle represents a filled and therefore unbuckled Si 
dimer. The triangles represent buckled dimers. The flat side 
of a buckled dimer indicates the upper Si atom with a filled 
dangling bond, whereas the pointed side indicates the lower Si 
atom with the empty dangling bond. The charge transfer from 
the Sr ad-atom to one of the surrounding dimers is indicated 
by the arrow, the preferred adsorption site (see section IVTll in 
the neighboring valley by the open circle. 



are in line with previous calculations for isolated Ba on 
Si(001)i2£ We believe that both atoms behave in a similar 
fashion. Ba has been found mostly on sites A in the 
trenches, but also on sites D on top of the dimer rows. 
The main difference between Ba and Sr lies in the energy 
difference between the two metastable sites A and D. 
For Ba the difference is 0.88 eV— which is substantially 
larger than the 0.29 eV for Sr. 

The chemical binding can be well understood in an 
ionic picture as suggested by the chemical binding of the 
silicides. The two electrons of the Sr atom are donated 
into an unoccupied dangling bond of a Si dimer. Interest- 
ingly we find this electron pair to be localized at a single 
dimer. This is evident from the dimer buckling, which 
vanishes when both dangling bonds are occupied. 

The Sr atom experiences an additional electrostatic 
stabilization from the remaining three buckled dimers 
next to it. They are buckled such that the negative, and 
therefore raised, silicon atoms are located next to the Sr 
atom. The local configuration is shown in Fig. 01 This 
arrangement significantly affects the buckling of the two 
dimer rows adjacent to the Sr atom: (1) the buckling gets 
pinned and is therefore observed in STM images in the 
vicinity of a Sr atom whereas it is thermally averaged 
out on the bare surface; (2) the dimer buckling within 
one row is reversed as already pointed out by Wang et 
ali^ This becomes apparent by looking at the row left of 
the Sr atom in Fig. 0] 

From the static structure shown in Fig. 0] it is not ev- 
ident why both dimer rows contain a buckling reversal 
as observed in Fig. 1 b) of Yao et ali^i We attribute the 
experimental observation to a dynamical effect: One elec- 
tron pair can rapidly migrate from the filled dimer to one 
of the three buckled dimers next to the Sr ad-atom. If 
this fluxional motion occurs on a time scale faster than 
the time scale for a buckling reversal of the entire chain, 
the buckling will appear pinned in both dimer rows. 

In our super cell with an even number of dimers in a 
row, every reversal of the buckling must be compensated 
by a second one, thus artificially destabilizing the site 
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Sr coverage 


Figure number 


Energy /Sr [eV] 


l/oo 


Ela) 


just schematic 


l/oo 


Hb) 


just schematic 


1/6 




-1.15 


1/6 


Eld) 


-1.15 


1/4 




-1.10 


1/4 


Eff) 


-1.06 


1/4 


im bottom 


-0.30 


1/4 


El top 


0.44 


1/2 


El 


-0.92 


2/3 


[TTlleft 


-0.74 


1 


ITUIleft 


0.04 


1 


ITUI right 


0.08 


1 


im right 


-0.12 


4/3 


m 


0.28 



TABLE IV: Energies per Sr ad-atom relative to our refer- 
ence energies for all structures graphically represented in the 
figures unless not already tabulated elsewhere in the paper. 



A. This adds an uncertainty of up to 0.12 eV to all en- 
ergies for the isolated Sr. Even taking this uncertainty 
into account, isolated Sr ad- atoms do not form a ther- 
modynamically stable phase at any coverage as will be 
demonstrated in the next section. 
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VI. CHAIN STRUCTURES AT DILUTE 
COVERAGES 

The Sr- atoms on the surface tend to arrange in chains, 
as seen in the STM experiments.-?*^ Similar results have 
been obtained for Ba on Si(001)i2^ 

Our calculations predict random, single chain struc- 
tures up to a coverage of 1/6 ML. Between 1/6 ML and 
1/4 ML we find condensed, single and double chain struc- 
tures. Above 1/4 ML the multiple chain structures con- 
vert into disordered arrays of double vacancies as will be 
discussed below. We investigated chain structures with 
coverages of 1/16, 1/10, 1/8, 1/6, 1/4, 3/10, 1/3, 5/14 
and 6/16 ML. We find an energy gain of 0.3 to 0.4 eV 
per Sr atom when single chains are formed from isolated 
Sr ad-atoms on A-sites. 

The chain formation is driven by the electrostatic at- 
traction between the positively charged Sr-ions, located 
at A-sites, and the negatively charged dimers: An iso- 
lated Sr ion located at an A-site, with the lone pairs of 
the four neighboring dimers pointing towards it, donates 
its two valence electrons into a silicon dangling bond ad- 
jacent to a neighboring valley. Thus it offers a preferred 
binding site for a Sr atom in that valley, namely next 
to this filled dangling bond as seen in Fig. 0) This sec- 
ond atom in turn donates its electron pair to the dimer 
row which does not already contain a negatively charged 
dimer and all four surrounding dimers will rearrange to 



FIG. 5: Schematic representation of a set of chain structures. 
Isolated chains are shown in the top row. The ordered single 
chain structures at 1/6 ML are shown in the middle row. The 
bottom row shows the ordered structures at 1/4 ML. The 
left hand side shows straight (diagonal) chains running at an 
angle of 60 % to the dimer rows, while the right column shows 
zig-zag chains perpendicular to the dimer rows. See Fig.|l|for 
a description of structural elements. The surface unit cells 
are outhned. The energies are listed in Tab.HVI 



point their lone pairs towards the new ad-atom. This 
process continues to form chains of Sr atoms. The filled 
dimers are clearly identified by the absence of any buck- 
ling. 

Two nearest Sr ions of one chain are displaced parallel 
to the dimer row by one lattice constant in order to po- 
sition the filled dimer inbetwecn. Since the favorable Sr 
positions are staggered with respect to the dimers, the 
chains run at an angle of 30 degree relative to the direc- 
tion of a dimer or 60 degree with respect to the dimer 
rows. The small energy difference of 0.02 eV per Sr ion 
between diagonal and zig-zag chains at e.g. l/8th ML 
cannot be considered as a hard number due to the sys- 
tematic errors of DFT calculations. It does, however, 
indicate that the energy cost for changing the direction 
of such a chain is negligible so that these chains may 
meander on the surface. 

As for the isolated Sr ad-atom, an additional stabiliza- 
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tion occurs due to the dimcr buckling of the surrounding 
sihcon dimers. The negatively charged, raised part of an 
adjacent dimer is located next to the Sr-ion, stabilized 
by electrostatic and covalent interactions. Reversing the 
buckling of one of the dimers next to a Sr ion raises the 
energy by 0.38 eVi^ This ordering induces a freezing of 
the dimer buckling, which reaches far out into the clean 
silicon surface, as can be clearly seen in the STM images 
by Hu et alm^ 

At first sight one might think that there is a second 
preferred binding site in Fig. ^ right next to the ini- 
tial ad-atom on the A site below the open circle. This 
configuration is, however, only possible for a pair of Sr 
atoms. A further continuation of such a chain perpendic- 
ular to the dimer row will make it impossible to rearrange 
dimers in a way that only filled dangling bonds are ori- 
ented towards the ad-atoms. Such chains arc therefore 
destabilized with respect to diagonal or zig-zag ones. 

A favorable registry between two Sr-chains is obtained 
if the Sr atoms are either in contact or separated by an 
even number of vacant A sites along each valley of the 
Si-surface. This follows from a simple building princi- 
ple which is an extension of what is already known from 
isolated chains: 

• Each Sr atom is electrostatically stabilized by four 
negatively charged silicon atoms located next to it. 
Negative silicon atoms have two electrons in their 
dangling bond and are in a raised, sp'^-like bonding 
configuration. 

Violation of this rule raises the energy by 0.38 cV 
per empty dangling bond next to the Sr atom. At 
this level of abstraction we do not distinguish be- 
tween a negatively charged Si atom of a buckled 
and an unbuckled dimer. 

• There are no reversals of the dimer buckling in the 
Sr-free regions on the surface. A buckling reversal 
increases the energy by 0.06 eV. This is a conse- 
quence of the anti-correlated coupling of the dimer 
buckling within one row (see section lrv)l . 

When the chains approach the shortest possible dis- 
tance before they collapse into double chains, we ob- 
tain a partially ordered structure at 1/6 ML as shown 
in Fig.|Hl(c) and (d). It should be noted that our calcu- 
lations indicate that also the condensed chains at 1 /6 ML 
change their directions frequently, even though synchro- 
nized with the neighboring chains running in parallel. 

As the coverage increases, Sr atoms arrange themselves 
into double chains as shown in Fig. EJc) , resulting in a 
partially ordered surface structure at 1 /4 ML. The " 1 x 2" 
areas in Fig. 4 of Ojima et ali^ can be explained by dou- 
ble chains at a 1/4 ML. Their interpretation that buckled 
Ba dimers at 1/2 ML coverage are responsible for this 
"wavy structure" cannot be supported by our calcula- 
tions. 

If we continue this building principle beyond 1/4 ML 
we will obtain sequences of triple, quadruple, etc., chains 



of Sr atoms separated by double vacancies. However, the 
positions of the double vacancies of different valleys are 
then only weakly correlated. Double vacancies of neigh- 
boring valleys can arrange themselves almost arbitrar- 
ily except that they do not line up perpendicular to the 
dimer rows. This implies a new building principle of dou- 
ble vacancies that do not necessarily arrange in chains. 
Note that a multiple chain structure of Sr atoms can also 
be interpreted as a chain structure of double Sr vacan- 
cies. This building principle can already be observed in 
Fig. [3fe) and (f) for 1/4 ML, which show that double 
vacancies in neighboring valleys can assume three out of 
four relative positions. While the energy difference be- 
tween the structures in Fig. Ete) and (f) is as small as 
0.04 eV per Sr atom, it decreases further for the anal- 
ogous structural patterns for wider Sr-chains. For cov- 
erages close to one-half ML, when the concentration of 
double vacancies is dilute, we therefore expect a nearly 
random arrangement of double vacancies instead of mul- 
tiple chains. 

The reciprocal lattice vectors for the diagonal chain 
structures are (|,0) and (±^,i), where n is the peri- 
odicity of the the real space unit cell (compare Figs.EJc) 
and (e)) along the dimer row direction. It should be 
noted that there is some structural disorder due to fre- 
quent changes in the chain direction. For zig-zag chains 
the corresponding reciprocal lattice is spanned by the 
vectors {j,Q) and (0, ^). The actual diffraction pattern 
observed in experiment will contain a mixture of both 
reciprocal lattices. 

The density of states of the single chain structure at 
1/6 ML is shown in Fig. The states on the Sr atom, 
which appear in the valence band, can be attributed to 
the tails of the Si dangling bonds, which hybridize with 
the Sr-s orbital. We observe states in the Si band gap, 
which are assigned to the empty dangling bonds on the 
buckled dimers. These states actually form a single band 
that is separated from the valence and conduction bands. 
The fact that they appear as individual states is an arti- 
fact of our discrete sampling of the Brillouin zone. This 
band pins the Fermi level in the lower part of the band 
gap. This feature remains nearly unchanged in the den- 
sity of states of the double chain structure at 1/4 ML 
(Fig. (Tj). It disappears, however, with the absence of 
the half-occupied dimers at 1/2 ML as seen in Fig. |H1 
Thus the gap states remain approximately in their posi- 
tion as the coverage increases, but the density of states 
is reduced. Hence the Fermi-level will remain pinned in 
the lower part of the silicon band gap up to a cover- 
age of 1/2 ML. At this point the Fermi-level becomes 
unpinned. In the following we will see that for coverages 
above 1/2 ML states from the conduction band are pulled 
into the band gap of silicon, and pin the Fermi level in 
the upper part of the band gap. 

This finding explains the discontinuity of the band- 
bending as observed by Herrera-Gomez et al^ Their 
XPS studies show that the Fermi level shifts up by al- 
most 0.5 cV when the coverage is increased from below 
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FIG. 6: Layer-projected density of states of the single chain 
structure at 1/6 ML. Values are divided by the number of 
atoms in the layer. The vertical line at eV indicates the 
Fermi level, the arrow points at the characteristic surface 
band in the bandgap of silicon. 



FIG. 7: Layer-projected density of states of the double chain 
structure at 1/4 ML. Values are divided by the number of 
atoms in the layer. The vertical line at eV indicates the 
Fermi level, the arrow points at the characteristic band states 
in the bandgap of silicon. 



1/2 ML to above. However, it should be noted that such 
a discontinuity in Fermi-level pinning is not specific to 
detailed structures: also higher energy structures exhibit 
a similar behavior. 



VII. THE (2 X 1) RECONSTRUCTION AT 1/2 
ML 

At a coverage of 1 /2 ML all dangling bonds of the sur- 
face dimers are fully occupied (Fig. O. It can be consid- 
ered as the canonical Sr-covered Si surface. It is the only 
Sr-covered surface structure without states in the band 
gap of silicon. This structure is " isoelectronic" to an 
hydrogen terminated silicon surface and is therefore ex- 
pected to be comparably inert. The increased resistance 
to oxidation has already been observed experimentallyii^ 
A (2 X 1) reconstruction at 1/2 ML has already been 
reported by several authorsiiiii^iii^ 

This structural template has been proposed to be the 
basic building block of the interface between silicon and 
SrTiOg-S 



VIII. FROM 1/2 ML TO 1 ML 

For coverages between 1 /2 ML and 1 ML, we find struc- 
tures built up from three structural templates: 

1. the (2 X 1) reconstruction at 1/2 ML (Fig. 13 

2. the (2 X 1) reconstruction at 1 ML (Fig. UUl left) 

3. the (1 X 1) reconstruction at 1 ML (Fig. UHl right) 

When we increase the coverage above 1/2 ML, the ad- 
ditional atoms occupy the D-sites, since all A-sites are 
already occupied. When all A-sites and all D-sites are 
occupied, as shown in the left panel of Fig^l the cover- 
age is that of 1 ML. At this coverage, two electrons are 
transferred to each silicon atom instead of only one as 
in case of the 1/2 ML. These electrons can fill the dimer 
antibonding states and thus break up the dimer bond. 
When all dimer bonds are broken at 1 ML we obtain a 
(1x1) reconstructed silicon surface with a Sr-ion above 
the center of each square of silicon atoms fFig. 1101 right). 
This structure is, however, never realized in its pure form 
due to the large strain in the top layer. An indication 
for the strain is the difference between the spacing of 
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FIG. 8; Layer-projected density of states of the (2 x 1) struc- 
ture at 1/2 ML. Values are divided by the number of atoms in 
the layer. The vertical line at eV indicates the Fermi level. 




FIG. 9: The (2x1) reconstructed surface at 1/2 ML coverage. 
The energy is listed in Tab. HVI 
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FIG. 10: Building blocks for the surface structure with IML 
coverage having a (2 x 1) (left) and a (1 x 1) (right) recon- 
struction. The outlined circles represent single Si atoms in 
the first layer. The energies are listed in Tab. II Vl 




FIG. 11: The (3 x 1) structures for coverages 2/3 ML (left) 
and 1 ML (right). The energies are listed in Tab.Hvl 



Sr atoms in bulk SrSi (Fig. and that of this hypo- 
thetical surface structure. The former is larger by 25 %. 
Nevertheless, this pattern is found as a building block in 
a number of low-energy structures between coverages of 
0.5 ML and 1 ML. The (1x1) reconstruction at 1 ML 
is unfavorable by 0.04 eV per Sr atom compared to the 
corresponding (2 x 1) structure. 

By combining the (2 x 1) structure at 1/2 ML (Fig.El) 
and the (1 x 1) structure at 1 ML (right panel of Fig. lTUIl . 
a scries of structures with increasing coverage and period- 
icity can be formed: dimer rows are separated by stripes 
of the (1 X 1)/1 ML structure with increasing width. This 
leads to a series of (n x 1) structures at a coverage of 
ML. We have investigated these structures from n = 3 to 
n ~ 6. 

The first structure in the series is the 3x1 recon- 
structed surface at a coverage of 2/3 ML shown in the 
left panel of Fig. ^2 It consists of dimer rows separated 
by a stripe of two Sr atoms in the (1 x 1)/1 ML config- 
uration. According to our calculations this structure is 
present as a distinct phase between 0.5 ML and 1 ML. 

In order to form this structure the dimer row pat- 
tern needs to reconstruct. This process is facilitated by 
the additional electrons in the conduction band which 
weaken the dimer bonds: beyond a coverage of 1/2 ML 
Sr is likely to act as a catalyzer for dimer bond rearrange- 
ment. 

The density of states for this (3 x 1) structure (Fig. ll2|) 
exhibits a surface band that ranges from the continuum 
of the valence band of bulk silicon into its bandgap. The 
surface band is localized on the dangling bonds of the 
silicon dimers and the undimerized silicon atoms. The 
Fermi level is pinned between this band and the conduc- 
tion band of silicon (see also discussion in section IvUl . 

For coverages greater than 2/3 ML, that is n > 3, the 
structures can be interpreted as stripes with coverage of 
1 ML with a (1 X 1) reconstruction separated by dimer 
rows without Sr, as discussed above. As the stripes with 
(1 X 1)/1 ML increase in width, they build up strain, 
which can be released by forming dimer rows with Sr 
atoms on top (Fig.^]left). This corresponds to a transi- 
tion between both 1 ML structures shown in FigllOl Since 
the energy difference between these structural variants is 
smaller than 0.04 eV per Sr atom, it is likely that they do 
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FIG. 12: Layer projected density of states for (3 x 1) surface at 
2/3 ML coverage. Values are divided by the number of atoms 
in the layer. Note that the characteristic peak just below the 
Fermi level is not a single state, but the top of a surface band 
extending into the continuum of the valence band states. The 
vertical line at eV indicates the Fermi level. 

not form distinct phases but transform into each other in 
a fluxional fashion. 

At a coverage of 1 ML we find a series of structures 
built from the two structural templates in Fig. 1101 In 
our calculations the (3 x 1) (Fig. [TT|l and (4 x 1) re- 
constructions are most stable and, within the theoretical 
error bar, degenerate. There is a slow increase in energy 
towards the (5 x 1) and (6 x 1) reconstructions which are 
less than 0.11 eV per Sr atom higher in energy. 



IX. BEYOND 1 ML 

Beyond 1 ML additional Sr atoms deposit on top of 
the silicon atoms in the (1 x 1) stripes. The first com- 
mensurate structure is the (3x1) structure at a coverage 
of 4/3 ML shown in Fig. US) 

Fan et al. have observed a (3 x 1) reconstruction at 
a coverage of 1.3 ML, Bakhtizin et aliSiS have published 
STM images with a (3 x 1) periodicity for a coverage 
of 1.2 ML. There, alternating bright and dark stripes 
have been observed, which we attribute to the Sr atoms 
on top of the dimer rows (darker stripes) and those on 
top of the of the two-fold coordinated Si atoms (more 




FIG. 13: 4/3 ML coverage showing a 3x1 silicide layer. The 
energy is listed in Tab. HVI 



prominent stripes). 



X. THE (3 x 2) STRUCTURE AT 1/3 ML 

Diffraction and STM studies identify a (3 x 2) recon- 
struction at 1/3 ML7iVi12,19i22,2.3,24,27,51^ ^hc diffrac- 
tion studies (LEED, RHEED) did not distinguish be- 
tween the orientations parallel and perpendicular to the 
dimer rows as they average over multiple terraces. Most 
previous studies assumed that the 3x direction of the 
(3 X 2) surface unit cell is orientated parallel to the dimer 
rowsiSiSiSLS However, the STM images of Hu et alm^ 
(Fig. 5) and Ojima et ali^i (Fig. 4) clearly show, that the 
3x axis is orthogonal to the dimer row direction. This is 
particularly evident from the images showing the phase 
boundaries between the (3 x 2) reconstructed domains 
and chain structures. This observation implies that the 
dimer row pattern is disrupted. 

Our lowest energy structures for this coverage are vari- 
ants of the quadruple chain as described in section IVII 
We have been unable to determine a thermodynamically 
stable structure with a (3 x 2) diffraction pattern at 
1/3 ML. 

Hu et alM- suggested two (3 x 2) structures for 1/6 and 
1/3 ML. We simulated both of them and found them 
energetically unfavorable compared to the correspond- 
ing chain structures by more than 0.61 eV per Sr atom. 
Despite various other attempts, we failed to arrive at a 
thermodynamically stable surface structure with this re- 
construction. 

This failure may be attributed to our inability to scan 
the entire phase space. However, one should also con- 
sider the possibility of co-adsorption of other elements 
such as hydrogen or oxygen, which may help to tie up 
the dangling bonds created by disrupting the dimer row 
pattern. These effects have not been considered in the 
present study. 

For the sake of completeness we also studied the model 
for the 1/3 ML coverage with the 3x direction parallel 
to the dimer row. It was lower in energy than the struc- 
tures suggested by Hu et alm^ Nevertheless it turned out 
0.23 eV per Sr atom higher in energy than the quadruple 
chain structure. 
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FIG. 14: Silicide structures with a c(4 x 2) reconstruction 
at 1/4 ML: The 1/4 ML c(4 x 2) structure as proposed by 
McKee et alp- (top) and the lowest energy structure we find for 
that chemical composition (bottom). The large dark circles 
correspond to Sr, the smaller, outlined ones to Si in the silicide 
layer (only in the top figure) and the small black dots to 
Si atoms in the first full layer. In the lower graph the Si 
atoms of the silicide layer have formed dimers indicated by 
the triangles. The energies are listed in Tab.HVI 



XI. SRSI2 SURFACE LAYER 
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FIG. 15: The adsorption energy— per (1 x 1) unit cell as a 
function of Sr coverage. The open diamonds represent ther- 
modynamically accessible structures, the triangles correspond 
to metastable structures. 



decomposition of two terraces into a single terrace with 
an additional half- monolayer of silicon. 

We consider the difference in formation energy of 
0.80 eV per Sr atom, relative to our lowest energy struc- 
ture at this coverage, as too large for this structure to be 
physically relevant. 



So far, we have discussed structures on a stoichiometric 
Si-surface. The atomic model for the interface between 
Si and SrTiOs by McKee et al.^ which can be seen in the 
top panel of Fig. has inspired us to also investigate 
reconstructions with only 1/2 ML of silicon in the surface 
layer. Such a structure can in principle be formed by the 
migration of Si atoms or dimers from step edges onto the 
terraces. For a clean surface this process is clearly not 
energetically favorable. However, it cannot be excluded 
a priori that the presence of Sr stabilizes a surface with 
half a monolayer of silicon. 

We started from the structure proposed by McKee et 
al^ which consists of a quarter monolayer of Sr and half 
a monolayer of Si. In this structure the Sr atoms occupy 
every second A-site in the valleys, while the Si atoms de- 
posit on top of the center of four subsurface Si atoms. We 
restricted our search to the c(4 x 2) periodicity reported 
by McKee et al. 

We find a number of metastable structures, the most 
stable one (Fig.^^bottom) differs substantially from the 
proposal by McKee et al. The half monolayer of silicon 
on the surface combines into dimers, a behavior already 
known from Si ad-atoms on Si(001)i^ The Sr atoms oc- 
cupy positions in the center of four such dimers. 

The energy of this structure is, however, still higher by 
0.80 eV per Sr ad-atom than our lowest energy structure 
at this coverage, namely double chains of Sr atoms. The 
energy was evaluated relative to the average of a 4-layer 
and a 5-layer slab, representing a terrace. This energy 
thus describes the process of adsorption of Sr and the 



XII. PHASE DIAGRAM 

We now investigate which of the reported structures 
form at given experimental conditions. The thermody- 
namic stability is determined by the zero-Kelvin Gibbs 
free energy G{fi) = E — ^N, where E is the energy per Sr 
atom and N is the number of Sr atoms, fj, is the chemical 
potential of the Sr atom relative to our reference energy 
for Sr. The extrinsic quantities, such as energies G and 
E, as well as the atom numbers are measured per (1 x 1) 
unit cell of the silicon surface. In Fig. ^| we show the 
adsorption energy E per (1x1) surface unit cell. 

The thcrmodynamically stable phases are determined 
by connecting the points in Fig. llSl bv line segments and 
forming the lower envelope. Each line segment corre- 
sponds to the coexistence of two phases, denoted by a 
and b, at the ends of the line segment with energies Ea 
and Et and Na and Nt Sr atoms per (1x1) unit cell. For 
a coverage A'' per (1 x 1) miit cell between Na and N^, 
the energy is E{N) = Ea + w^^i^ ~ ^a)- The slope 

of the line segment, namely /x = '^^J^^ = ^"Z^l ' 
chemical potential at which both phases a and b coex- 
ist. All structures not contributing to the lower envelope 
are not thcrmodynamically stable at low temperatures. 
At higher temperatures they may be stabilized due to 
entropic effects. 

The one-dimensional phase diagram is shown in Fig llBI 
The region in between two lines correspond to the differ- 
ent surface phases. The lines indicate the coexistence of 
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FIG. 16: One-dimensional phase diagram of Sr on Si(OOl) as 
a function of the Sr chemical potential. The shaded region 
between the 1/4 and 1/2 ML phase indicates disordered sur- 
face structures as described in section IVII Above the dashed 
line at a chemical potential of eV bulk Sr silicides are ther- 
modynamically stable compared to surface structures. 



to occasional Sr vacancies, which help to release strain. 
It seems surprising that no phase boundary is seen at 
1/2 ML. This may be due to the fact that on the one 
hand the (2x1) structure develops continuously out of 
the multiple chains structures below 1/2 ML and on the 
other hand it can be transformed continuously into the 
(3x1) structure by introducing thin stripes with local 
coverage of 1 ML. 

The multiplicity of structures with low energy above a 
coverage of 2/3-ML suggests the presence of disordered 
structures at elevated temperatures. McKee reports in- 
commensurate structures beyond a coverage of 5/8 MLii& 

A similar multiplicity of structures is found for 1 ML. 
Our calculation predicts a (3 x 1) structure as the most 
favorable. However we find also (4 x 1), (5 x 1) and (6 x 1) 
reconstructions within a window of 0.11 eV per Sr atom. 

We did not extend our calculations beyond 4/3 ML, 
where an overlayer of metallic Sr is formed. Therefore, 
our data does not necessarily indicate the presence of a 
phase boundary at 4/3 ML. 



the two neighboring phases. 

We find thermodynamically stable phases at 1/6 ML 
and 1/4 ML. We did not fully explore the phase dia- 
gram below 1/6 ML. In this region we find single chain 
structures as shown in Fig[51(a) and (b). Entropic effects 
disorder the arrangement of chains at room temperature. 

At 1/4 ML we find an ordered structure of double 
chains. Beyond 1/4 ML we predict nearly random ar- 
rangements of double vacancies. The double vacancy is 
stabilized relative to the single vacancy due to a favorable 
dimer buckling. 

The structure at 1 /2 ML is a clear phase boundary for 
a wide range of chemical potentials. 

The next phase boundary is found at a coverage of 2/3 
ML with a (3 X 1) reconstructed surface. Above 2/3 ML 
there are a number of low-energy structures with various 
coverages and a (n x 1) periodicity. 

Within the phase region of the 2/3 ML coverage, bulk 
Sr silicides become thermodynamically stable as indi- 
cated by the dashed line in Fig.^] As mentioned above, 
we expect the onset of silicide formation to be signif- 
icantly delayed due to thin-film effects as the bulk Sr 
silicides are highly incommensurate with the Si(OOl) sub- 
strate. 

It is of interest to compare our phase diagram with 
that obtained by McKee et ali^ At around 600 C they 
determine three line-compounds at 1/6 ML, 1/4 ML and 
5/8 ML. Our calculations reproduce phase boundaries 
at 1/6 ML and 1/4 ML. The next phase with a (3 x 1) 
diffraction pattern is seen at 5/8 = 0.625 ML, which is 
close to the coverage of 2/3, where our calculations pre- 
dict a phase with the identical periodicity. The differ- 
ence in coverage corresponds to a change by one in 16 
ad-atoms or 1/24 ML. This difference may be attributed 



XIII. CONCLUSIONS 

In this paper we have investigated the surface struc- 
tures of Sr adsorbed on Si(OOl) as a function of coverage. 
We propose a theoretical phase diagram by relating the 
phase boundaries at zero temperature to chemical poten- 
tials, which can be converted into partial pressure and 
temperature in thermal equilibrium. We predict phases 
at 1/6 ML, 1/4 ML, 1/2 ML, 2/3 ML and 1 ML. Struc- 
tural models are discussed for all experimentally observed 
reconstructions except a (3 x 2) reconstructed layer at- 
tributed to a coverage of 1/3 ML. The models are ex- 
plained in terms of structural templates and rationalized 
in terms of their electronic structure. 

Our findings elucidate the chemistry of alkaline earth 
metals on Si(OOl) and the phases of Sr on Si(OOl), which 
is expected to provide critical information for the growth 
of one of the most promising high-K gate oxide to date, 
namely SrTiOa. 
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